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A new type of Langevin equation exhibiting a non trivial 
phase transition associated with the presence of multiplicative 
noise is introduced. The equation is derived as a mesoscopic 
representation of the microscopic annealed Ising model (AIM) 
proposed by Thorpe and Beeman, and reproduces perfectly its 
basic phenomenology. The AIM exhibits a non-trivial behav- 
ior as the temperature is increased, in particular it presents a 
disorder-to-order phase transition at low temperatures, and a 
order-to-disorder transition at higher temperatures. This be- 
havior resembles that of some Langevin equations with multi- 
plicative noise, which exhibit also two analogous phase tran- 
sitions as the noise-amplitude is increased. By comparing the 
standard models for noise-induced transitions with our new 
Langevin equation we elucidate that the mechanisms control- 
ling the disorder-to-order transitions in both of them are es- 
sentially different, even though for both of them the presence 
of multiplicative noise is a key ingredient. PACS: 05.40.+j 



I. INTRODUCTION 

A great deal of attention has been recently devoted to 
the study of physical effects induced by the presence of 
noise, i.e. phenomena appearing in stochastic systems, 
which would be absent in the sole presence of the de- 
terministic part of the corresponding Langevin equation 
PJ. By now it is clear that noise can generate quite un- 
expected and counterintuitive behaviors as, for example, 
the stochastic resonance in which the output to in- 
put ratio of a bistable system subjected to the presence of 
an oscillating force is strongly enhanced by the presence 
of an additional stochastic term. Other examples are 
the resonant activation ||] , and the noise induced spatial 
patterns Q. Another type of phenomena the noise is at 
the base of, are the so called noise induced phase tran- 
sitions. These came to light in an interesting paper by 
Van der Broek, Parrondo and Toral [[| (see also Jl|,^,0|). 
These authors pointed out the fact that some Langevin 
equations may exhibit a noise- induced ordering transition 
(NIOT), i.e. a phase transition that is not expected from 
the analysis of the deterministic part of such equation. 
The phenomenology is as follows: 

i) For low enough noise amplitudes the system is dis- 
ordered (i.e. the order parameter takes a zero value). 

ii) At a certain critical value of the noise amplitude the 
system exhibits a NIOT and, in a range of noise intensi- 
ties above it the system remains ordered. 



iii) Finally, for noise amplitudes larger than a second 
critical value, the noise operates in a more standard way, 
this is, disordering again the system. We refer to this sec- 
ond phase transition as noise induced disordering transi- 
tion (NIDT). 

A physical explanation of the NIOT was proposed in 
H ; The ordering of the system is the consequence of the 
interplay between the noise and the spatial coupling 
In particular, the noise generates a short time instability 
at every single site, and the presence of a spatial coupling 
renders stable the non trivial state generated in that way. 

A minimal model capturing the essence of the NIOT 
has been recently proposed |J . It has been clarified that 
the essence of the NIOT is purely multiplicative, this is, 
in order to generate an ordering transition the noise has 
to appear multiplied by the field variable. In this way, 
it has been possible to recognize that the NIOT is char- 
acterized by a set of critical exponents other than those 
of well established universality classes (as for example, 
that of the Ising model) |u|!|. It has also been shown 
that due to the multiplicative origin of this transition it 
is possible to observe the phenomena in d = 1, dimension 
at which is very unusual to observe phase transitions. 

Other results concerning NIOTs and NIDTs can be 
found in the literature [pd|-p^[. A common feature of all 
the previously referred models, is that they are defined by 
means of Langevin equations, this is, equations describ- 
ing the physics at a mesoscopic, coarse-grained scale (in 
fact, the concept of noise is meaningful only at this level). 
In this context, it is an interesting task that of analyzing 
microscopic models that exhibit similar non-trivial be- 
haviors; one of which is the anneal Ising model J^J. By 
studying the connection between a microscopic systems 
and their respective mesoscopic representation one could 
shed some light on the way in which microscopic mecha- 
nisms generate the very non-trivial effects described at a 
mesoscopic scale. 

In what follows we introduce the time honored anneal 
Ising model. It was proposed and described more than 
twenty years ago by Thorpe and Beeman |l3j ]. A more 
detailed description of it will be presented in the next sec- 
tion; here we summarize the main properties we are in- 
terested in. The system is an Ising model in which the in- 
teractions, J, among spins, are annealed (not quenched) 
random variables that change from bond to bond and 
are extracted from a fixed probability distribution, P(J). 
Under certain conditions (this is, for some distributions 
P(J) to be specified later), the system phenomenology is 
as follows: 
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i) For low temperatures the system is disordered, i.e. 
the averaged magnetization is zero. 

ii) At a critical value of the temperature, T-y, the system 
exhibits a second order phase transition. As the temper- 
ature is further increased above T\ the averaged mag- 
netization keeps on growing until it reaches a maximum 
value and it starts decreasing if T is increased further. 

hi) At a second critical temperature, T%, the system 
exhibits another phase transition (analogous to the well 
known ferromagnetic-paramagnetic disordering transi- 
tion of the standard pure Ising model). The system re- 
mains disordered for temperatures higher than T2. 

This phase diagram resembles very much the behav- 
ior of the previously described noise induced transitions 
in Langevin equations. It is our purpose here to analyti- 
cally derive a coarse-grained, mesoscopic, representation, 
in terms of a Langevin equation, of the microscopic an- 
nealed Ising model (AIM) to further explore the eventual 
relations between both phenomena. 



II. THE ANNEALED ISING MODEL 

Let us consider a d-dimensional impure Ising model in 
the sense that the value of the coupling constant among 
spins, J, changes from bond to bond, being an annealed 
random variable with a fixed temperature-independent 
probability distribution, P(J) (which is not quenched but 
annealed at every site) . Following the strategy proposed 
by Thorpe and Beeman jl4| the model can be exactly 
mapped into a standard pure Ising model with an effec- 
tive parameter, K — J/T, that depends on P(J) and T, 
and we write as K e f f (T) . In particular E3] , 



dJ 



P(J) 



coth[K eff - J/T] - e(K eff ) 







(1) 



where t{K) is the correlation function of two nearest- 
neighbor spins in the pure Ising model. By solving 
the implicit equation Eq. (Q) one obtains K e ff as a 
function of the temperature and the parameters char- 
acterizing P(J). Note that, in particular, for the two- 
dimensional case, the Onsager's solution Jl5[ provides an 
explicit value for e(K) and therefore Eq. ([j])~can be solved 
and, furthermore, the system magnetization can be ex- 
pressed as a function of T. Let us suppose now that, in 
particular, the distribution P{J) is centered at a posi- 
tive value Jo (favoring ferromagnetic ordering), and has 
a variable width (standard deviation), 5 J. The resulting 
magnetization for this particular type of distribution is 
qualitatively represented in figure 1 (see also BM). 




FIG. 1. Magnetization as a function of the temperature T, 
for different values of S J; = SJi < 8J2 < SJs < 8J4 for the 
annealed Ising model. 

Observe that for narrow distributions of J the magne- 
tization curve is similar to its corresponding in the pure 
Ising model. Instead, as 6 J is increased, a disordering 
tendency is observed at low temperatures, and in partic- 
ular, for a large values of the width (as for example, 5 J4 in 
Fig.l) the system is disordered at low temperatures, and 
exhibits a disorder-to-order phase transition at a certain 
temperature. The standard ferromagnetic-paramagnetic 
(order-to-disorder) transition is also present and occurs 
at a variable value of T for different values of 8 J. 

The physical mechanism leading to the previous be- 
havior was argued in Q to be the competition between 
ferromagnetic and antifcrromagnetic types of interactions 
that emerges when sufficiently large values of 6 J are con- 
sidered. In particular, when 6 J > Jo both positive and 
negative values of the coupling constant are accessible 
at each bond, and in that case , for low temperatures 
the system is in a frustrated state in which ferromagnetic 
and antiferromagnetic domains compete. That frustra- 
tion makes the ferromagnetic order parameter to vanish. 
As the temperature is further raised the thermal noise 
activates annihilation of domain walls and the system is 
more likely to ordinate. As a result, the averaged magne- 
tization grows with increasing temperature. At a given 
point this effect ceases, and the standard role of the tem- 
perature as a disorganizing source sets to work. 



III. CONTINUOUS REPRESENTATION 

Let us now follow a standard procedure |l(| to cast 
the previous AIM into a continuous Langevin equation. 
For that purpose we first consider the pure Ising model 
case, and write down its associated equilibrium partition 
function: 



Z 




(2) 



Introducing auxiliary Gaussian integrals in terms of con- 
tinuous variables <pi (with i varying from 1 to the total 



2 



number of spins, N, in the lattice), and performing the 
change of variables ipi = K^ifij we obtain (l6) 



Z oc / dipi...di(jisrexp 



i 

(3) 



Expanding the hyperbolic-cosine in power series, per- 
forming a transformation to Fourier space, considering 
only the leading dependence on the temperature, and 
performing the continuous limit we finally obtain |16| 



Z oc d[ijj] e 



-H 



H = ljd d x [K Q (1 - 2K )^ 2 (x) + P(*K - l)(V^f 



(4) 



with K Q = J d d xK{x) 1 and p=l/2j d d xK(x)x 2 . In this 
way we have derived a Ginzsburg-Landau coarse grained 
Hamiltonian for the Ising model. This could have been 
guessed a priori by using heuristic arguments, but we 
have preferred to follow the previous procedure that per- 
mits to obtain explicit expressions for the coefficients as 
function of the microscopic parameters. In this way, ob- 
serve, for example, that both the diffusion constant and 
the coefficient of the quadratic term depend on the cou- 
pling through Kq, therefore in order to simplify the no- 
tation we define the diffusion constant, D = p(4Ko — 1). 
Taking only the main relevant dependences on D we can 
write 



H 



d d x 



aD 9 b . D 



(5) 



where a is a tuning parameter proportional to the dis- 
tance to the critical temperature, and 6 is a positive pa- 
rameter. Let us stress once more that we are neglecting 
higher order dependences of b and a on D, and we assume 
them to be unessential to reproduce the microscopic phe- 
nomenology of interest at mesoscopic level (this hypoth- 
esis will be verified afterwards). The simplest Langevin 
equation with a stationary distribution characterized by 
a Gibbsian distribution with the Hamiltonian in Eq. (|5|) 
is well known to be |l7| , [l8|] 



dtip = -{aD + bijj 2 )i> + DV 2 ip + n(t) 



(6) 



where r){t) is a Gaussian white noise with (r](x,t)) = 0, 
and (77(2, t)rj(x',t')) = 5 d (x - x')S(t - t'). 

At this point we can analyze the effects of an annealed 
distribution of J in the microscopic AIM at the level of 
Langevin equations. For that purpose let us observe that 
in order to mimic the variability of the coupling in the 
AIM we can just substitute D at each site in Eq. (||) 
by a stochastic variable, namely: D — > D + £(x,i), with 
(£(x,t)) = and (Z(x,t)Z(x',t')) = a 2 D S d (x-x')S(t-t'), 



where D and <td play the role of Jo and S J respectively 
in the microscopic model. In this way we obtain, 

d t ip = -[a(£> + + 6V7 2 ]^ + 7JV 2 i/. + V(CV^) +»?(*) . 

(7) 

This equation (intended in the Ito interpretation [[18]) 
constitutes our continuous representation of the AIM. Let 
us underline that the differences with respect to the pure 
case, Eq. (|^, are two: the presence of a multiplicative 
noise, and an extra term that couples spatial fluctuations 
of D with V"0. Changes of a, parameter which appears 
multiplying both the linear term and the multiplicative 
noise, correspond to temperature variations. 

We have analyzed Eq. (]?]) in mean field approximation 
|p|,^9|Js|], and by performing systematic numerical simu- 
lations in two dimensions. The mean field approxima- 
tion is performed along the lines discussed in ||^,p^fl . 
For the numerical simulation wc have employed the Eu- 
ler method p9[ , in a 32*32 lattice, with lattice spacing 
Aa = 1, and considered a time mesh At = 0.001. With- 
out lost of generality the parameters b and D have been 
fixed to 1 and 10 respectively. Different noise amplitudes, 
ctd, have been considered. The main results we have ob- 
tained are as follows: in both, mean field approximation 
and in the numerical simulation, we reproduce the qual- 
itative behavior of the order parameter as a function of 
the temperature characteristic of the microscopic model 
(see Fig. 2 and Fig. 3 and compare them with Fig.l). 
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FIG. 2. Averaged magnetization of the Langevin equation 
as a function of a (i.e. the temperature), in mean field ap- 
proximation, for different values of an- 

In mean field approximation the order-to-disorder crit- 
ical point is located at a = 0, and in numerical we obtain 
also a close to zero critical value which does not depend 
on a/). On the contrary, the location of the disorder- 
to-order transition, depends on 073, analogously as the 
location of T\ depends on 6 J in the AIM. Observe that 
this transition is not sharp in the lowermost curve of 
Fig. 3 due to finite size effects. Curves in Fig. 2 and Fig. 3 
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change with increasing o~d in the same way they do in 
the AIM when increasing SJ, i.e., the larger the noise the 
smaller the ordering. 




FIG. 3. Averaged magnetization ol the Langevin equation 
as a function of a (i.e. the temperature), in numerical simu- 
lations, for different values of o\d. 

Let us stress once more that in order to obtain the 
transition, we change both the coefficient of the linear 
term and of the multiplicative noise term. If one of 
these two coefficients was kept fixed while the other was 
changed the microscopic phenomenology would not be re- 
produced. The presence of the multiplicative noise term 
is essential to generate the disorder-to-order transition. 
We have performed a numerical study of Eq. (Jt|) omitting 
the term proportional to V(£V^>) , and conclude that 
none of the previous conclusions is qualitatively affected 
by this suppression; by omitting this term the disorder- 
to-order critical point is shifted to a lower value of a, and 
consequently this term has only a disorganizing effect. 
We could consequently write down a minimal model just 
by dropping out this unnecessary term, in the same way 
we omitted other irrelevant higher order dependences on 
D in the derivation of the Langevin equation. We con- 
clude that the proposed Langevin equation with multi- 
plicative noise in the Ito representation reproduces qual- 
itatively all the interesting properties of the anneal Ising 
model, and in particular the reentrant phase transition. 
Therefore, once more it is shown that the multiplicative 
noise is the key ingredient of highly non-trivial phenom- 
ena in stochastic systems at a mesoscopic level. 

Let us finally remark that the phenomenon we have 
just described is not the usual noise-induced transition 
as reported in previous works First of all, in 

those works only the multiplicative noise amplitude has 
to be changed to obtain a NIOT, while in our case the 
transition is obtained by varying the parameter a that 
multiplies both the multiplicative noise and the linear 
term. Consequently in our case the disorder-to-order 
transition is not purely noise-induced. Second, con- 



sidering a Stratonovich representation of the Langevin 
equation with multiplicative noise is essential in those 
works to generate noise-induced ordering. In fact, stan- 
dard Langevin equations as those described in p,0,|| 
do not exhibit NIOTs when intended in the Ito repre- 
sentation ( [^|). On the other hand, in the model pre- 
sented here, the Langevin equation is intended in the Ito 
sense, and due to its peculiar structure, namely the cou- 
pling between a and that we have justified from a 
microscopic point of view, it can exhibit a rather rich 
phenomenology. In particular the system shows an or- 
dering and a disordering transition as the temperature 
is increased but it does not exhibit, for example, the 
short time instability characteristic of the phenomena 
discussed in I 
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